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We develop a boundary field theory approach to electron spin resonance in open S = 1/2 Heisen- 
berg antiferromagnetic chains with an effective staggered field. In terms of the sine Gordon effective 
field theory with boundaries, we point out the existence of boundary bound states of elementary ex- 
citations, and modification of the selection rules at the boundary. We argue that several "unknown 
modes" found in electron spin resonance experiments on KCuGaFe [I. Umegaki et ai, Phys. Rev. 
B 79, 184401 (2009)] and Cu-PM [S. A. Zvyagin et ai, Phys. Rev. Lett. 93, 027201 (2004)] can be 
understood as boundary resonances introduced by these effects. 
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Introduction. - Impurities often introduce new as- 
pects in physics, Kondo effect being a notable example. 
In particular, impurity effects in strongly correlated sys- 
tems are currently among central topics in condensed 
matter physics. Although the standard perturbation the- 
ory can fail, there arc number of powerful theoretical ap- 
proaches to strongly correlated systems, especially in one 
dimension. Impurity effects in gaplcss one-dimensional 
systems have been vigorously studied in terms of bound- 
ary conformal field theory. In contrast, impurity effects 
in gapped one-dimensional systems received much less 
attention, with the exception of the edge states in the 
S = 1 Haldane gap phase [1-3]. 

On the other hand, integrable models and field the- 
ories have been successfully applied to many gapped 
one-dimensional systems. In quantum magnetism, a 
field-induced gap in 3 = 1/2 Heisenberg antiferromag- 
netic (HAFM) chains is described in terms of a quan- 
tum sine Gordon field theory [4-6]; the one-dimensional 
Ising chain with critical transverse field and a weak lon- 
gitudinal field realizes a quantum field theory with Eg 
symmetry [7, 8]. Experimental studies indeed found el- 
ementary excitations predicted by these integrable field 
theories. Application of integrable field theories to impu- 
rity/boundary effects in gapped one-dimensional strongly 
correlated systems is an interesting but largely unex- 
plored subject. 

In this Letter, we present a theory of electron spin 
resonance (ESR) in S = 1/2 HAFM chains in a stag- 
gered field with boundaries, which may be realized by 
nonmagnetic impurities. The low-energy effective the- 
ory of the system is the quantum sine Gordon field 
theory with boundaries. In fact, this theory is inte- 
grable even in the presence of a boundary, and bound- 
ary bound states (BBS) of elementary excitations have 
been found in the exact solution [9-14]. The existence 
of BBS, and modification of the selection rules, imply 
extra resonances in addition to those in the bulk. ESR 
measurements on corresponding systems KCuGaFe [15] 
and [PM-Cu(N0 3 )2(H20)2]n (PM denotes pyrimidine, 



abbreviated as Cu-PM) [16, 17] had found several res- 
onances that could not be accounted for by the theory. 
We argue that several of those resonances can be suc- 
cessfully identified in terms of the boundary sine Gordon 
field theory. 

Boundary sine- Gordon field theory. — We consider a 
semiopen chain with the Hamiltonian 



n = [ JS i ■ S ^ -^bH- S -S j + {-l)°D ■ Sj x Sj-i] , 

j<0 



(1) 

which models one side of an infinite chain broken by a 
nonmagnetic impurity at j = 0. g is the g tensor of 
localized spins, and /ib is Bohr magneton. 

The Zeeman energy can be represented as 
-MiiE^^li + (-1)^4 Hereafter, we 
assume that g^ b gS a b and |<7^ b | <C g, and employ a unit 
h = k B = gfJ-B = 1. 

We consider H = Hz applied along the z direction (z 
is a unit vector in the z direction). The last term of (1) 
is the staggered Dzyaloshinskii-Moriya (DM) interaction. 
This can be eliminated by a staggered rotation of spin 
about the direction of D by angle (— lp'a/2, where a = 
tan _1 (|D|/J) where j is the site index. 

Under an applied field, this transformation leaves a 
staggered field h ~ H x D/2, which is perpendicular to 
H and to D. [5, 6] Together with the staggered field due 
to the staggered component of the g tensor, the effective 
model may be given by 



H 



j<0 



S, ; - US] - h(-l) j Sf 



(2) 



keeping only the most important terms, h = c s H is the 
effective staggered field, approximately perpendicular to 
the applied field; the direction of the staggered field is 
chosen to be the x axis. c s depends both on the staggered 
DM interaction and on the staggered part of g tensor g s ab . 



Using bosonization formulas, 
1 



S ~. 



2irR 



d x (j) + C*(-l) x cos(<j)/R + Hx), (3) 



n ±2-rrRi 



*[C^(-lT + C^co S (<t>/R + Hx)}, (4) 



at low temperature T <C J, the model (2) is mapped to 
a boundary sine Gordon (BSG) field theory, defined by 
the action 



A 



dt 



dx 



C^/icos (2ttR4>) 



(5) 



v is the spin-wave velocity. The fields <f> and <fi are dual 
and compactified as <j> ~ + 27ri? and ~ <p + 1/R; 
m is the uniform magnetization density and the nonuni- 
versal constants s and C^ s arc determined numeri- 
cally [18]. The bulk operator cos(27ri?0), which repre- 
sents the transverse staggered magnetization, is relevant 
in the renormalization group sense. Thus it induces a 
finite mass (excitation gap), as it was observed in the 
experiments. [5, 6] The bulk gap is stable against dilute 
nonmagnetic impurities. 

In the absence of the staggered field term, the field 
4> obeys the Dirichlet boundary condition <f>(x = 0,t) = 
const. [19]. This is equivalent to the Neumann boundary 
condition in terms of the dual field: 



d x (j>(x 7 t) 



0. 



2=0 



(6) 



Inclusion of the staggered field could change the bound- 
ary condition; in fact, the staggered field in the bulk 
would also induce the corresponding operator cos (2TrR(f>) 
at the boundary. If this boundary perturbation is dom- 
inant, <p would obey the Dirichlet boundary condition. 
However, since cos (27ri?0) is, as a boundary operator, 
(nearly) marginal in RG, its effects are negligible at the 
energy scale set by the bulk spin gap. Thus, for a small 
staggered field h, the boundary condition can still be re- 
garded as the Neumann on </>. 

Energy spectrum. — The elementary excitation of the 
bulk sine Gordon field theory includes soliton (denoted by 
S) and antisoliton (S) with the same mass M. A soliton 
and antisoliton carry soliton charge Q = +1 and Q = — 1, 
respectively. Additional particles called breathers are 
generated as bound states of a soliton and an anti- 
soliton [20, 21]. There can be several different kinds 
of breathers B n with the mass M„ = 2M sin(n7r£/2), 
for n = 1, 2, . . . , [£ -1 J • Here |_ - J is the floor function. 
Breathers have zero soliton charge. The soliton charge 
is conserved in the bulk sine Gordon field theory. For 
the case of our interest, that is experimental situations in 
KCuGaFg and Cu-PM, the parameter £ = 1/(2/ttR 2 - 1) 
satisfies £ < 1/3. In particular, ^ 1/3 in the low field 
limit H 0. 



Fortunately, the BSG theory (5) with the Neumann 
boundary condition (6) is still intcgrablc [9]. An anal- 
ysis of the boundary S matrices implies [9-12, 22] the 
existence of the BBS with the mass 



A/bbs = M sin(7r£), 



(7) 



for £ < 1/2. Therefore, the BBS with A/ B bs ~ y/SM/2 
does exist in the low field limit of the present system. 
The BBS of soliton, antisoliton, and first breather turn 
out to be identical and there is only one type of BBS. 
This is another manifestation of the soliton charge non- 
conservation at the boundary. 

Thus the analysis of the BSG theory predicts a new 
excited state, which is a BBS, at the energy Mbbs lower 
than the bulk gap M . In fact, in an earlier numerical 
study of an open chain based on density matrix renor- 
malization group, Lou et al. [23, 24] had found such 
an excited state localized near the boundary. They 
called it a midgap state. Figure 1 shows a compar- 
ison of the soliton (12) and BBS (7) masses, to the 
numerically obtained bulk gap and the energy of the 
midgap state in Refs. [23, 24]. Here we used M sa 
1.85(/i/J) 2 / 3 [ln(J//i)] 1 /6 appropriate for H = 0, which 
was used in Refs. [23, 24]. The excellent agreement 
between the two energies means that the midgap state 
found in the DMRG calculation was nothing but a BBS. 

While it was pointed out in Refs. [23, 24] that the 
midgap state is localized near the boundary, physical un- 
derstanding of its origin has been lacking. The analogy 
to the edge state in the Haldane phase is clearly inap- 
propriate, since the 5 = 1/2 HAFM in a staggered field 
is topologically trivial. With the present identification of 
the midgap state with the BBS, the BSG theory is es- 
tablished as an effective theory describing the boundary 
physics of the system. This is essential in understanding 
the ESR spectra. 

Electron spin resonance. — Next we discuss the ESR 
spectrum, which is given by I(u>) oc ^Xphysi^ = Oi^- 1 )- 
Here, x" is an imaginary part of the dynamical suscep- 
tibility x- Xh — (q^) = ~ Gs+s- (9, ioj = oj + ie) where 
positive infinitesimal e is the analytic continuation of the 
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Figure 1. Numerically obtained bulk gap (open circles) and 
midgap (solid circles) in [23] are compared with the soliton 
mass M and the BBS mass A/bbs (7). 
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temperature Green's function Gg+s- {Q> The stag- 
gered rotation of spins by angles (— iya/2 to eliminate 
the DM interaction mix the uniform (q = 0) and the 
staggered (q — n) components. The physical susceptibil- 
ity Xp hys (q = 0,u) is 



(8) 



where D z is the z component of the DM vector D parallel 
to the applied field. On the right-hand side, we dropped 
the longitudinal susceptibility Xzz(l = ^,^) because it 
contains the same resonances as those in the transverse 

part Xh (<Z = 0, uj), and merely modifies their intensities. 

Let us first review the ESR in the bulk, in the limit of 
T — > 0. The uniform part x+-(0,w) reflects transitions 
caused by the operator 



J q=0 



^cos (<P/R + Hx). 



(9) 



The operator cos (4>/R) changes the soliton charge by ±1, 
and thus must create at least one soliton or antisoliton. 
The other factor e ±27rfl - 1 * can create any number of exci- 
tations with zero soliton charge in total. Thus ESR in- 
duced by the operator (9) with the lowest possible energy 
corresponds to creation of a single soliton or antisoliton. 
It should be also noted that the factor Hx in the cosine 
causes the shift of the momentum with H. That is, the 
created soliton or antisoliton should carry the momentum 
H . Thus the ESR due to a single soliton/antisoliton cre- 
ation is at frequency uj = E$ = v 7 M 2 + H 2 . There are 
also resonances due to (9) at higher energies, correspond- 
ing to creation of additional elementary excitations. 

Next we turn to the staggered part X+-(l = f^), 
which reflects transitions caused by the operator 



9* 

^ q— 71 



(10) 



This carries zero soliton charge, and thus the simplest 
excitation induced by this operator is creation of a single 
breather B n . This leads to the resonances at w = M n . 

Now let us discuss the boundary effects on ESR. Here 
we discuss the ESR spectrum based on physical picture, 
leaving systematic formulations to the Supplementary 
Material . 

First we consider the contribution of the staggered part 
X+-(<7 = tt,lo). The simplest excitation created by the 
operator (10) is a single breather. In the presence of the 
boundary, the first breather can form the BBS. Thus the 
resonance with the lowest energy in the presence of the 
boundary is given by U) = Mbbs ■ Creation of a breather 
not bounded at the boundary and the BBS is also possi- 
ble, leading to the resonance at u) = Mqbs + M n . 

In order to understand the boundary effects on ESR, 
we need to clarify the issue of the momentum conserva- 
tion. In general, total momentum is conserved due to 





Bulk 


Boundary 


T — 


U = E S ,E S +Mn 


oj = E n , E s + Mbbs, E n + Mbbs 


T > 


u = \E S - M„\ 


w = E n — Mbbs 



Table I. Typical resonance modes in x+-(l = 0,uj). Reso- 
nances shown in the second row are absent at T = 0. The 
soliton resonance Es is accompanied with all bulk resonances. 
On the other hand, -Es does not necessarily appear in the 
boundary resonances. In fact, some boundary resonances are 
involved with a novel resonance E n instead of Es- 





Bulk 


Boundary 


T = 


uj = M„,M„ +M m 


uj = Mbbs, M„ + Mbbs 


T > 


uj = M n - M m (n > to) 


uj = M n - A/bbs 



Table II. Typical resonance modes in x+-(Q — ~k,uj). 
Breather masses are directly measured in the bulk part. The 
BBS mass itself appears in the boundary resonances. Equa- 
tions (8) shows that intensities of these modes are smaller 
than those in Table I by (D z /J) 2 . 



the translation invariancc of the system. In the pres- 
ence of the boundary, the translation invariance is lost 
and the total momentum is no longer conserved. Nev- 
ertheless, the momentum is still important in discussing 
ESR spectra, because the momentum of each elementary 
excitation is conserved or reversed in a scattering with 
another elementary excitation, or in a reflection at the 
boundary. Thus, once created, the set of momenta of el- 
ementary excitations is conserved, up to the sign of each 
momentum. 

The existence of the boundary has an interesting ef- 
fect, in addition to the contribution of the BBS, on the 
ESR spectrum. Although the expectation of the operator 
cos (<j)/R) vanishes in the bulk, it is nonvanishing [25, 26] 
near a boundary with the Dirichlet boundary condition 
on 4> [equivalent to the Neumann boundary condition (6) 
on <fr]. This reflects <\> taking a fixed value at the bound- 
ary. Thus, in the vicinity of the boundary, the leading 
contribution from the uniform part is effectively given 



by the operator S ~ e 



' cos (Hx) , which creates 



excitations with zero soliton charge, in contrast to the 
original one (9). This is another consequence of violation 
of soliton charge conservation at the boundary. The cre- 
ated excitations should carry the total momentum ±H , 
up to the uncertainty ~ l/l p . Here l p is the pinning 
lengths cale, namely <fi(x) ~ const, if |a;| < l p . 

Thus, the simplest among the possible ESR processes 
due to this operator is the creation of a single breather B n 
with momentum ±H. This corresponds to the frequency 



Mr, 



H 2 . 



(11) 



The resonance at E n with n > 1 was absent in the bulk 
and is a new feature due to the boundary. We emphasize 
that, these new resonances do not simply follow from 
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Figure 2. (color online) Frequency vs field diagrams of ESR in KCuGaFe [15] for (a) H \\ a and (b) H || c configurations. The 
dotted line is high temperature paramagnetic resonance ui = H. Open symbols and filled symbols represent bulk modes and 
"unknown modes." Es,Bi,B2,B3 denotes resonances lj = Es , Mi, M2, M3 by a soliton S and breathers B n . An antisoliton S 
also leads to ui = Es- The labels Ui, U2, U3 are "unknown" peaks found in [15]. (a) The configuration H || a bears the smallest 
h — c s H with the coefficient c s — 0.031. An excitation u> = Mi + A/bbs is found in addition to the bulk excitations, (b) We 
have the largest staggered field, c B = 0.178, when H |] c. This large c s makes the rich kinds of boundary modes detectable. 
The labels BB, 2*Bi-BB, B2-B1+BB denote to = Mbbs, 2Mi - M B bs,M 2 - Mi + M B bs- 



the existence of the midgap state numerically found in 
Refs. [23, 24]. In fact, the resonance frequency E n does 
not explicitly contain the energy Mbbs- This shows the 
necessity of the BSG framework to fully understand the 
physics at the boundary. 

At finite temperatures, additional resonances may be 
observable. When the initial state contains the BBS as 
a thermal excitation, a resonance at uj = M n — Mbbs 
exists, corresponding to annihilation of the BBS and cre- 
ation of a breather B n . Similarly, when the initial state 
contains B\ , the resonance at w = M n — M\ + Mbbs cor- 
responds to the creation of B n and binding of B\ at the 
boundary. These resonances are contained in the stag- 
gered part x+-{' k t lo )- The uniform part also contains, 
at finite temperatures, additional resonances. 

Typical resonance modes are summarized in Tables I 
and II. Note that intensities of resonances due to the 
staggered part x+-{q = and x'Lil = n ^ w ) are SU P~ 
pressed by the factor (D z /J) 2 in (8), compared to those 
from the uniform part x'+-(q = 0,7r). 

Comparison with experiments. — Thus several novel 
resonances, which are absent in the bulk, are derived from 
the BSG theory. They indeed match the "unknown" reso- 
nances observed previously in ESR spectra on KCuGaFg 
(Figure 2) and Cu-PM (Figure 3). Here the uniform 
field effect is taken into account in the soliton-antisoliton 
mass [27], 



M 



2v T(g/2) 

v^r((i + o/2) 



r(i/(i + £))c s J 



(l+«)/2 



(12) 

While the ratio c s = h/H can in principle be determined 
by the staggered DM interaction and g tensor, we use the 
value obtained by fitting experimental data. 

Each figure shows different resonances, though. The 



difference between the spectra in two materials can be 
understood in terms of the different magnitude of the DM 
interaction. The coefficient c s at its maximum is larger 
(0.178) in KCuGaF 6 compared to 0.083 in Cu-PM. While 
precise estimates of the DM interaction and the staggered 
g tensor are not available, the staggered DM interaction 
is presumably larger in KCuGaFg. This leads to larger 
mixing of the staggered part X+— (9 = 7r i a -')- Thus it is 
natural that the resonances due to the mixings are ob- 
served only in KCuGaFg (Fig. 2). We note that, the 
simplest possible BBS contribution at u = Mbbs is not 
observed for H || a in Fig. 2 (a). This is presumably 
because only two of the frequencies used in the experi- 
ments can detect the resonance below the bulk resonance 
at bj = Mi. On the other hand, in the Cu-PM with a 
smaller DM interaction, only the contributions from the 




5 10 15 20 
Magnetic field (T) 

Figure 3. (color online) Frequency vs field diagrams of ESR in 
Cu-PM [16]. Two "unknown" peaks U\ and U2 are attributed 
to uj — E2, Ei + Mi — Mbbs respectively, where E n is defined 
in Eq.(ll). These boundary modes appear in XxxiQ = 0, to). 
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uniform part are observed (Fig. 3). We conjecture that, 
with more careful examination of the spectra, more res- 
onances due to the BBS will be found in experiments. 

Conclusions. — We point out the existence of BBS and 
modification of the selection rules at boundaries of the 
5 = 1/2 antifcrromagnetic chain in an effective staggered 
field, in terms of the BSG theory (5). The boundary 
effects can account for the mysterious "unknown modes" 
found in two compounds KCuGaF 6 [15] and Cu-PM [16, 
17]. 

In the compound KCuGaFg, magnetic ions Cu 2+ and 
nonmagnetic ions Ga 3+ form a pyrochlorc lattice. Mag- 
netically, the compound KCuGaFg is effectively regarded 
as S = 1/2 HAFM chains of Cu 2+ ions. However, since 
Cu 2+ and Ga 3+ ions occupy equivalent positions, an in- 
tersite mixing of them can occur in the course of synthe- 
ses [28]. We speculate that the intersite mixing brings 
about nonmagnetic impurities in spin chains. We ex- 
pect that a few percent of the nonmagnetic impurities 
would lead to observation of the boundary resonances as 
discussed in this Letter. Our picture may be verified ex- 
perimentally by controlling the density of nonmagnetic 
impurity, which will change the intensity of boundary 
resonances. 
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Supplemental Material to 
"Boundary Resonances in S — 1/2 Antiferromagnetic Chains under a Staggered Field" 



EFFECTIVE LATTICE MODEL 
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Figure 4. Semi-infinite Heisenberg antiferromagnetic chain 
with an effective staggered fieid, cut by a non-magnetic im- 
purity (the gray filied circie). 

We derive the effective lattice model (2) of the main 
text, with an emphasis on the mixing of q = tt modes 
through the transformation eliminating the DM interac- 
tion. First we consider the model Hamiltonian 



E 

j<0 



[JSj ■ Sj-i - HS~ + (-l) j D ■ Sj x Sj-i] . (13) 



Here, for simplicity, the effects of the anisotropic g tensor 
are ignored, and they will be taken into account later. We 
rotate the spin Sj about the direction of D by the angle 
(-lp'a/2 where a = tan _1 (|D|/ J). When |U| < J, 
a <C 1 and the transformation of the spin operator can 
be linearized as 



^ + (-1)^x5,, 



(14) 



Thus, in the presence of the uniform applied field H , a 
staggered field perpendicular to both H and D is gener- 
ated. 

In addition, the staggered component of g tensor also 
contributes to generating the staggered field. In contrast 
to the one due to the staggered DM interaction, the stag- 
gered field coming from the staggered component of the 
g tensor is not necessarily perpendicular to the applied 
field. However, in the effective field theory description, 
the longitudinal component of the staggered field accom- 
panies an oscillating factor, which makes it vanish upon 
spatial integration. Thus it can be ignored for the pur- 
pose of discussion of gap opening, and only the transverse 
staggered field has to be considered. 

Therefore, the S — 1/2 antiferromagnetic chain under 
an applied magnetic field, with both staggered g tensor 
and staggered DM interaction can be described by the 
effective lattice model (2) of the main text. The effec- 
tive staggered field h depends on both the staggered DM 
interaction and the staggered g tensor. 

In the course of eliminating the staggered DM inter- 
action, the spin operator is transformed according to 



Eq. (14). This has to be taken into account in the dis- 
cussion of physical susceptibility. Namely, the physical 
uniform susceptibility Xph y s(9 = 0j w ) nas contributions 
from both uniform and staggered susceptibilities calcu- 
lated for the effective model after the elimination of the 
DM interaction: 



XphysO? = 0,w) ~ x+-(q = 0,w) 

' A x 2 
J 



X+-(<? = 7T,W) + 



J 



x'Lii = 



(15) 



where D± is the magnitude of the component D per- 
pendicular to the z axis. As we will see below, the lon- 
gitudinal susceptibility x'Lii = 7r 7 w ) adds no resonant 
peak to the former two susceptibilities. The longitudinal 
staggered susceptibility just modifies intensities of the 
transverse susceptibilities. 

Fig. 4 shows a schematic picture of the staggered fields 
generated from the external uniform field H, in the bulk 
(h) and at the boundary (/is). Note that the boundary 
staggered field hs, coupled to the end spin Sq, is indis- 
tinguishable to the bulk staggered field h in the lattice 
model. But their behaviors under renormalization group 
(RG) transformations are quite different. 



RG ANALYSIS 

The bosonization of spins leads to an effective bound- 
ary field theory with an action, 



1 f 1 



A= dt dx 

J — oo J — oc 

- C^hcos(2nR0) 



- C^h B / dt cos [2ti^ R(f>(x = 0)] . 

(16) 



Scaling dimensions of the bulk and boundary staggered 
magnetizations are respectively 



Xh = ttR , 



2irR z 



(17) 



When the uniform field H is zero, the effective lattice 
model (2) in the paper is SU(2) symmetric, where the 
compactification radius R is given by R = 1/v2tt. As 
the uniform and staggered fields increase, the compact- 
ification radius decreases [1]. Recall that a bulk opera- 
tor is relevant (irrelevant) when its scaling dimension is 
smaller (bigger) than 2, and a boundary operator is rel- 
evant (irrelevant) when its scaling dimension is smaller 
(bigger) than 1. Thus, both the bulk (xh < 1/2) and 
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hg = +00 




R 



Figure 5. The boundary RG flow. On the line hs = de- 
noted by "N", the boundary condition is Neumann. On the 
line hs = +00 denoted by "D", the boundary condition is 
Dirichlet. 




Figure 6. (a) A reflection process at the boundary. The 
boundary is represented as an infinitely heavy particle B 
(thick gray line), (b) A formation of BBS using a breather 
Bi. The system is excited from the ground state |0) = £>|vac) 
to the one-BBS state |1). 



boundary {xh B < 1) staggered magnetizations are rele- 
vant. The boundary staggered magnetization is marginal 
at zero field H — 0. The boundary RG flow is shown in 
Fig. 5. The line Kb = corresponds to the Neumann 
boundary condition. In the presence of the boundary 
staggered field, the boundary condition is neither Neu- 
mann nor Dirichlet: 



d x if>(x,t) 



2-KRC^h B sin[27ri?(^(a; = 0,t)] (18) 



liB = +00 corresponds to the Dirichlet boundary con- 
dition because the field <j>(x = 0,t) is pinned to a cer- 
tain value in order to optimize the potential energy 
calls cos[2nR(f>(x = 0, t)]. 

As we perform the RG transformation iteratively, the 
bulk staggered field h grows rapidly. When the effec- 
tive coupling h/J becomes 0(1), the RG transformation 
breaks down. On the other hand, the boundary staggered 
field Kb I J is still much smaller than 1 after the itera- 
tive RG transformations because the boundary staggered 
magnetization is almost marginal. Therefore, at the low- 
est order approximation, we may ignore the boundary 
staggered field h. Finally we reach the effective action of 
the boundary sine Gordon model (5) together with the 
boundary condition (6), in the main text. 



BOUNDARY BOUND STATE 

Here, following Rcfs. [2, 3], we briefly review the energy 
spectrum of the boundary sine Gordon model with the 
Neumann condition. The energy spectrum in bulk is the 
same as those in the conventional sine Gordon model. A 
soliton (denoted as S) and an antisoliton (denoted as S) 
exist and several breathers are formed as bound states 
of one soliton and one antisoliton. One can know that 
masses of breathers by analyzing simple poles of the exact 
S matrix. As is well known, the S matrix characterizes 



two-particle scatterings: 

ZaMZaM) = s££(*i - e 2 )z b2 {6 2 )z bl {(>i) 
zl{0i)zl{8 2 ) = sS&o?! - e^zl^zlie,) 

Z ai {Q 1 )Zl 2 {e 2 ) = 2ir5 aia2 8{6 1 - B 2 ) 

+ s b a i a b i(e 1 -e 2 )zl(9 2 )z bx (e 1 ) 

Operators Z a {6) and Z\{d) are called as Faddeev- 
Zamolodchikov (FZ) operators. The parameter is called 
as a rapidity, which parameterize the energy E a and the 
momentum P a of a particle with an index a as 

E a {9) = M a cosh 0, P a (6) = M a sinh 6. 

Since every scattering in intcgrablc field theories is fac- 
torizable to several two-particle scatterings, the S ma- 
trix determines whole energy spectrum in bulk. An n-th 
breather has a mass, 



M n = 2M sin 



(19) 



for n = 1, 2, • • • , [C _1 J ■ Since breathers are formed by a 
soliton and an antisoliton with the degenerate mass M, 
the mass of the breather M n is smaller than 2M. 

Similar procedure is applicable to investigate the en- 
ergy spectrum at the boundary. The scattering at the 
boundary is shown in Fig. 6 (a). The spatial boundary 
is regarded as an infinitely heavy particle B. 



Z\(Q)B = RU6)ZI(-6)B 



(20) 



The index a represents the anti-particle of the particle 
a. The reflection factor R\{0) has been studied in de- 
tail [3] . Simple poles of the reflection factor R h & (8) repre- 
sents energy levels localized at the boundary, called as a 
boundary bound state (BBS). 

For later convenience, we consider a rotated reflection 
factor 



K ab {9) 



^R b J^- 



ITT 

2 



(21) 
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According to Rcf. [3], the rotated reflection factor K ab {9) 
has the following simple poles 9 = iu a b within the physi- 
cal strip < Im0 < 7r/2, 



u a b = u S s = tt£, (a, b = S,S) 



(22) 



7T 7f£ 

■"ab = u Bi b x = 2 - "2"' (a = 6 = Si). (23) 

Here we assumed the experimentally realized situation 
3 < C" 1 < 4. Both poles (22) and (23) lead to the same 
BBS. In fact, the excitation gap of the BBS, Mbbs, is 
derived by the following two manners. 

1. One soliton S or one antisoliton S forms a BBS: 

Mbbs = M sin u ss = M sin(7rf ) . (24) 

2. One first breather B\ forms a BBS: 

M B bs = Mi smu Bl B 1 = Msin(7r£). (25) 

In (25), we used the relation (19). The formation of 
a BBS is sketched in Fig. 6 (b). Note that the BBS 
is formed by a soliton with the mass M or by a first 
breather with the mass Mi. The mass of the BBS should 
be smaller than min{M, Mi}. Since min{M, Mi} is the 
lowest excitation gap, the BBS appears below the bulk 
gap. 



where E and P are the total energy and the total momen- 
tum respectively. The state \B) is called as a boundary 
state written in FZ operators and the rotated reflection 
factor [2]. 



\B) 



e K lvac) 



exp 



JET ^K ab (0)zi(-e)zl(9) 



vac) 



(27) 



The Fourier transform of the x-ordcrcd correlation 
function is given by 



drdX e^ aT - q ^G s+S -{xux 2 ,T). 



(28) 



where Q is the Matsubara frequency. Instead of spatial 
variables x\ and x%, we used r = — \x\ — x%\ and X = 
(xi +X2)/2. The infinitely large N represents the length 
of the semi- infinite spin chain and the coefficient 1 /N can 
be understood as a density of non-magnetic impurities. 
When we compute dynamical susceptibility, we perform 
the analytic continuation iQ — > uj + ie, where e is positive 
infinitesimal. 



CORRELATION FUNCTIONS 

We summarize correlation functions in the presence of 
the boundary. In the Euclidean space, a two-point corre- 
lation function Gg+g- (x%, x%, t) depends on two spatial 
variables x\ and X2 because the translational symmetry 
is broken by the boundary (impurity). For a systematic 
treatment, it is useful to make a Wick rotation to regard 
the spatial variable x as the "time" , and t as the "space" 
variable. In this framework, 



G 



s+s 



-(xi,x 2 ,t) 



(v ! xc\T x S+(x 1 ,t)S-(x 2 ,0)\B) 
(vac|B) 



(26) 



The state |vac) is the ground state (with length given 
by the inverse temperature of the original system, and 
with the periodic boundary conditions). The symbol T x 
denots the x ordering, that is, the ordering which orders 
operators with larger x^, to right. Note that the spin 
5 a (x,r) is related to 5" (0,0) via 



S a (x,r) =e- x ^e~ lT ^S a (0,0)e 



JtP xP 
& 1 



BULK RESONANCES 



Here we derive several resonant peaks in the ESR spec- 
trum. An n,-particlc state \6\, 62, • • ■ , 9 n ) aia2 ... ar and its 
conjugate an "' a2ai (9 n , 



, 62, 9i I are defined by 



llMa,- •• ,e n ) aia2 ... an =Zt i (9 1 )Zl(9 2 )---Zl(9 n )\v a .c), 
Q "- Q2Q1 (0„,... ,0 2 ,0i| = (vacl^J^)---^^)^^). 

FZ operators Z\{9) and Z a {9) respectively creates and 
annihilates a particle with a rapidity 9 with an index a = 
S, S, B m . The energy E = M a cosh 9 and the momentum 
p = M a sinh 9 are parameterized by the single parameter 
9. We put v = 1 for simplicity. The above mult i- particle 
states are orthonormal. 



a ™'" ai (6»„, • ■ • ,9i\9[,--- ,Oai-o^ 
= S nm S ai< ■■■5 an< (27r) n 5{9 1 -6[)-- 

A set of n-particle states is complete: 



(29) 



I vac) (vac I 



1 y 



V 

^ n. 

ri—l a\ 



d6\ • • • d0 n . 
(2tt)« ' 



; w n/ a\ ■ --an 



(30) 
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This identity operator allows one to expand correlation functions by intermediate states. We can expand the boundary 
state in the power of K as 



e^lvac) 



I vac) + K\vac) 



(31) 



Let us consider the bulk part of the correlation Gs+s-(q — 7t,zcj), that is, 



G ( °l s _(q = TT,iui) = ^-C^ dr dr dX e lQT (vac\T x e^ R ^ Xl ^e^ R ^ X2 ' 0) \va,c) 



N 

1 VI m 



1 



1 



iui — M m iui + M n 



(32) 



Thanks to the Lorentz invariancc of the sine Gordon 
model, the matrix element (yac\e~ l27rR ^ Xl ' T ^\9) a involv- 
ing the intermediate state \9) a can be computed easily. 



(vac e 



-i27rR(f)(xi ,t) I 



-i2^H0 iT M m sinh 6»+xiM m cosh < 
r B m e 



taincd, which are respectively 



M 



D± 
J 



2 nvCf 



M 



{0 O )r6(w-Eg), 

cl>/R + i2-KR4> 



s 



(0 )\ J 6(w-Es). 



(35) 
(36) 



The factor FJ.J " = (vac|e- i27ri? ^ '°) \6) Bn is a form 
factor of the operator e,- %2vR v . Form factors generally 
depend on the rapidity 9, but the F^^^ 4 " does not. 
In integrable field theories, form factors can be exactly 
derived. One can find several exact form factors of the 
sine Gordon model in Ref. 4. 

After taking analytic continuation iui — > ui + ie, we 
obtain 



ImG 



(0) 

s+s- 



. (q = 7r, uj + ie) 



= E 

n=l 



\F\ 



?***\ 2 8(!*-M n ) 



(33) 



Multi-particle resonant peaks such as 5(ui — M n — M m ) 
exist in the ignored terms. (33) leads to peaks in (15) as 



6» = In [{H + VM 2 + H 2 } /M] . Both of them detect the 
soliton resonance at ui = E$ = VKP + H 2 . The latter 
(36) has a weaker intensity by the factor (D±/J) 2 than 
the former (35). The q = u component of S*, 

S*~C*(-l) x cos(cl>/R + Hx), 
is similar to the q = component of S^r, 



Thus, all resonant peaks in x'z Z (q = appear also in 
X+-(q = Oj w )i an d the latter has stronger intensities in 
(15). We may forget x'Lil = n i UJ ) when we focus on 
resonance frequencies 

We should emphasize that the breather resonances 
S(u! — M n ) can result only from the transverse staggered 
component x+_(<7 = in Xp hys (9 = 0,w). This is 

an important consequence of the mixing (15) caused by 
the staggered DM interaction. Without the mixing, we 
cannot observe any breather resonance. 



J 



n=i 



ft 



'\ 2 S(uj~M n ) + 



(34) 



BOUNDARY RESONANCES 

Transverse staggered component G s + s -(q = n,iui) 



As stated in the paper, the intensity of the staggered 
susceptibility x'+ - (<Z = t, is suppressed by the fac- 
tor (D z /J) . Bulk parts of the resonant peaks which 
comes from x"_(0,oj) and x Z2 (7r,w) arc similarly ob- 



Boundary scatterings induce rich amounts of addi- 
tional resonances. For instance, we consider the first or- 
der contribution in the expansion (31) to the G$+ s- (q = 

it, iui), 
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i poo pu i poo in 

= ±Cy dr drdXe^-J" — K a \0){v&c\T x e- l2 * R ^ xl ^e l2 * R ^^\ - 9,9) ab (37) 

N JO J-oo 2 7b"'" oo27r 

JV JO J-oo ^ 0jf> J-oo ^ ITT 

x [^(x 2 -a;i)(vac|e- i2ffii ^ 3!l ' r )|^)B 1 ' Bl (^|e i27r ^ a:2 ' )|-^e) a6 
+ H (xi - x 2 ){vac\e i2 ^ x ^\e 1 ) Bl Bl Q^e*****^ \ - 9,9) ab ] + ■ • ■ (38) 

9h(x) is the Heaviside step function. As we reviewed above, the factor K ah (9) has simple poles in the physical strip 
< Im0 < tt/2, 



Res K ab {9) = -if s , (a, b = S, S), (39) 

8=iuss 

Res K ab {9) = -if Bl , (a = b = B 1 ), (40) 

e=iu Bl b 1 

where u$s and u Bi b 1 are (22) and (23). The residues —ifs and —ifs x are n °t derived at the moment of publication. 
Here we assume fs,fB ± > 0. In the expansion of (37), the most dominant intermediate state is one-breather state 
|0i)b„. The singularity of K ab {9) affects (37) when and only when n—1 and a = b = B\. 

i 9 r°° 

G« g _ (q = tt, iw) ~ -C^ j dr j drdX e^/ Bl * | 2 

sin iabjBj cos ubiBi i pi~Mi sin ub 1 b 1 ^2XMi cos ub^! \ _i_ 

_^ 2 / Sl |^™IV i i 



2A/i cos UBxBi \iu) — MisiniiBiB! iui + Misuiub 1 b 1 
The pole Mi sin ub 1 b 1 = M sin(7r£) = A/bbs gives a resonance of the boundary bound state (BBS) in Xph ys (l = 0; w ) : 

— ) A^sin^/2) 6{UJ ~ MBBS) (41) 

The intensity of the peak (41) has two features, (i) The intensity is suppressed by the factor (D±/J) 2 . (ii) The 
intensity is independent of a concentration (1/N) of non-magnetic impurities. The first feature is common with the 
other resonant peaks which originates in x+-(q = tt,u>). The second is also common with the other bulk resonant 
peaks. However, resonant peaks involved with the BBS generally depends on the concentration. For instance, let us 
consider a two-particle intermediate state \9i,02)b 1 b 1 in Ggl s _(q = n,iuj) instead of the one-particle state |#i)bi in 
(38). The corresponding term in Gg+ S _ (q = tt, ilj) is 



rfr drdXe^- K B ^(9)- 



N s Jo J-oo 2 < /_ 00 2tt ^'2^ (2tt)2 

Here we assumed £i < £2 for simplicity. This process leads to a resonant peak in Xphys(9 = 0,w) as follows. 

' ' s J 1 BlBl 1 1 U V — -S(oj - Mi - Mbbs) (42) 



N \ J J 2M 1 2 sin(7r£/2) 

The two-particle form factor FjQf*(0i ~ #2) = (vac|e i2,rfl *( o ' o )|0i, 6> 2 ) BlBl is positive when 6»i - 9 2 = 
—iuB 1 B 1 i ~iuB lBl + *7T [4]. Note that the intensity is proportional to the concentration of non-magnetic impuri- 
ties because of the factor 1/N. 

The resonance frequency to — Mi + Mbbs is & simple superposition of two resonances oj = Mi and u> — A/bbs- If 
(2) 

we consider G s + S _ (q = tt, iuj), we will obtain multi-BBS resonances such as u> = 2MbbSj M n + 2Mbbs- 
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Transverse uniform component G s + s -(g = 0, ioj) 



Below we point out that some boundary resonances cannot be understood by any superposition of the existing 
resonance frequencies. This new resonance appears in x"_(y = 0,w). When we are concerned with q as component 
of correlation function Gg+s- (ffjiw), we may replace 

S£ ~ ~Y [Ose~ iHx + O s e lHx ] . (43) 

The fields O s = e - i( - 27rR $+<f>/R) and 0\ = e i(?**$+<1>/R) can create n solitons, m antisolitons and I breathers with 
n—m = 1 and n—m = —1 respectively Here I can be arbitrary. Thus, minimal intermediate states in G5+ 5- (q = 0, iuj) 
are a one-soliton state |#i)s and a one-antisoliton state \0i)§- 

The uniform component \+ - (q = 0, uj) does not lead to the resonance at uj = Mbbs- This is because a "shift of the 
wavenumber" is accompanied with x+-(q — 0,w). The bulk soliton resonance occurs at u> = E$ = \JM 2 + H 2 , not 
at uj = M . Es is obtained by substituting q = ±H into the Lorentz invariant dispersion relation E{q) = \J M 2 + q 2 
of the soliton. On the other hand, the BBS does not have y-dependent dispersion. Thus, the resonance u = A/bbs 
cannot originate in x"_(<z = 0, ui). 

An interesting resonance appears in processes involved with multi-particle intermediate states. Let us consider a 
case in which a two-particle state \6\, 62) b±s appears as an intermediate state. 

1 c 2 r r° ^v^iv- r r 



■ dr drdXe l " T -Y / — K ab {6) \ 



2 

1h{x2 - xi)e- iHr (va.c\O s (x 1 ,T)\9 1 ,6 2 ) BlS SB i (6 2 , fli|0^(x 2> 0)| - 6, 9) ab 

+ H (x 1 ~x 2 )e- iHr (v^\Os(x 2 ,0)\6 1 ,e 2 ) BlS SBl (e 2 ,9 1 \Ol(x 1 ,T)\ -9,9) ab 



Both B\ and S in the intermediate state \9i,6 2 )b 1 S can 
form the BBS. When the first breather B\ forms a BBS, 
the resultant resonance frequency is 

uj = M cosh #0 + M sin U55 . 

Since, M cosh6» = VM 2 + H 2 = E s and M sinu S s = 
A/bbS; this resonance frequency is equal to 

lo = E s + Mbbs- (44) 

This is a simple superposition of the bulk resonance u> = 
E$ and the boundary resonance uj = A/bbs On the other 
hand, when the soliton S forms a BBS, the resultant 
resonance frequency is somewhat strange, which is 

uj = E!+M B bs- (45) 
The first term is given by 

E 1 = ^M 2 + H 2 . 

This resonance frequency (45) cannot be represented 
by any superposition of existing resonance frequencies. 
Since a soliton S is trapped at the boundary and trans- 
formed into a BBS, a first breather B\ propagates in the 
bulk as if it were a soliton. 



INTENSITIES OF "UNKNOWN PEAKS" 

In our paper, we identified the "unknown modes" 
found in recent experiments on KCuGaFg [5] and Cu- 
PM [6, 7] with the boundary resonances. Here we briefly 
comment with intensities of the "unknown peaks" . The 
compound KCuGaFg has larger staggered field (c s = 0.18 
for H || c) than Cu-PM (c s = 0.083). Thus, in Cu-PM, 
intensities of the staggered susceptibilities x'+- (<7 = uj) 
and Xzzil = n ^ U! ) are strongly suppressed because the 
strength of the staggered DM interaction is equal to 
c s « 1 at most. We may expect that, if boundary 
resonances exist, they are likely to come from the uni- 
form susceptibility = 0,cj). In fact, our assign- 
ments of the "unknown modes" in Cu-PM, which are 
uj = Eb 2 ,Eb 1 + All ^ A/bbS; arc consistent with this 
speculation. One can find that the breather resonances 
uj = M n are weaker than the "unknown mode" U± in 
Ref. 6. We assigned U-y to u = Eb 2 - While it looks coun- 
terintuitive that the boundary mode has the stronger ab- 
sorption intensities than those of bulk modes, this rela- 
tion can be understood by the weak mixing of the stag- 
gered susceptibilities in (15). The breather resonances 
come from the staggered part. 

KCuGaFg has relatively stronger staggered field. 
Ref. 5 shows that the breather resonances at uj = M n 
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have approximately the same intensities as the soliton 
resonances uj = Eg. This experimental observation is 
also consistent with (15) and the large c s . As for the 
"unknown modes" , our assignments are consistent with 
the estimation of the large c s . We assigned "unknown 
modes" to w = M BBS ,Mi + M BBS ,2Mi - M BBS ,M 2 - 
M\ + MbbSj ah of which originate in the transverse part 
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